Random Reverse Cyclic matrices and screened harmonic oscillator 



Shashi C. L. Srivastava^Q and Sudhir R. Jain 



•0 



CO 

o 

CN 

43 
<D 
P-h 

CN 



43 ■ 

43 : 



> 

CN 
CN 

o 

CO 



1 RIB Group, Variable Energy Cyclotron Centre, 1/AF Bidhan nagar, Kolkata 700 064, India 
2 Nuclear Physics Division, Bhabha Atomic Research Centre, Mumbai 400 085, India 

(Dated: February 13, 2013) 

We have calculated the joint probability distribution function for random reverse cyclic matrices 
and shown that it is related to an iV-body exactly solvable model. We refer to this well-known model 
potential as a screened harmonic oscillator. The connection enables us to obtain all the correlations 
among the particle positions moving in a screened harmonic potential. The density of nontrivial 
eigenvalues of this ensemble is found to be of the Wigner form and admits a hole at the origin, in 
contrast to the semicircle law of the Gaussian orthogonal ensemble of random matrices. The spacing 
distributions assume different forms ranging from Gaussian-like to Wigner. 

PACS numbers: 05.40.-a, 05.45.Mt, 03.65.Ge 



I. INTRODUCTION 

Connections between random matrix theory and ex- 
actly solvable models are very important and interest- 
ing [jJ-Q- It is well-known that the invariant random 
matrix ensembles are related to some exactly solvable 
many-body problems in one dimension, as was found by 
Calogero and Sutherland In particular, the joint 

probability distribution function (JPDF) of the eigenval- 
ues of random matrices shares the functional form with 
the probability density corresponding to the quantum 
ground state of the TV-body problem. This observation 
is important as it allows one to obtain the correlation 
functions of one problem by knowing those for the other, 
comparing terms using a dictionary. In the same vein, 
even for the explanation of intermediate statistics Q, a 
random matrix model was found (Io| which, in turn, was 
related to an iV-particle system with an inverse-square, 
repulsive two-body interaction, and, an inverse-square, 
attractive three-body interaction 0, Q • Even for pseudo- 
Hcrmitian Hamiltonians (where there exists a metric rj 
such that H* = rjHr] -1 ), a random matrix theory can 
be built [ill; Ell- The connection of this with exactly 
solvable models is explored in Q. In turn, the models 
found in and @,i 

can be mapped to integrable 
[lH and chaotic systems [14[ for which, quite remarkably, 
analytically exact eigenfunctions are obtained. 

In this work, we study random, reverse cyclic matrices, 
that are real symmetric, 
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I ai a 2 ... djv^ 
a 2 a 3 ... at 



(1) 



/ 



with matrix elements chosen from an appropriate distri- 
bution function. Bose et al. [l5| derived the limiting 
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spectral distribution for reverse-cyclic matrices, but the 
JPDF and the spacing distribution function remain open 
problems. In fact it will be interesting to see how the 
special symmetric matrices having a very small degree 
of freedom (only N in this case), differ from the results 
known for their counterparts having a full degree of free- 
dom [i.e. N ^ N 2 +1 ^ ]. We have recently obtained the JPDF 
for the cyclic matrices, which forms another example in 
which the degree of freedom of the matrix is constrained 
(again only N) in an asymmetric matrix for which a spec- 
trum of a spacing distribution from Gaussian, to Wigner 
and not-so- Wigner type |16J is obtained. The results were 
also used to study a random walk problem on a one- 
dimensional disordered lattice [l?} where the evolution 
matrix is cyclic. On the one hand, there is a vast liter- 
ature about different results for random cyclic matrices 
in literature (see [l|| [l9| etc.). However, the same is not 
true for random reverse-cyclic matrices. Interestingly, 
reverse cyclic matrices appear (albeit with the name re- 
verse circulant and retro-circulant) in models for particle 
masses, flavour mixing, and CP violation. Here fami- 
lies of particles can be shown to emerge by a sponta- 
neous breakdown of discrete Z@ chiral symmetry, by the 
Higgs sector [2(J. The presence of reverse-cyclic matri- 
ces is due to S3 cyclic permutation symmetry of the La- 
grangian. Quoting Adler, "...in the limit of S3 cyclic per- 
mutation symmetry, we shall find that the fermion mass 
matrices in both the three and six doublet models are 
retrocirculants..." [2fjj. In another instance, while explor- 
ing whether discrete flavor symmetry S3 can explain the 
pattern of neutrino masses and mixings, reverse-cyclic 
matrices (again referred to as retro circulant) have been 
used as a perturbation matrix[2lj. It was also shown in 
[2lj that after third order perturbation, neutrino mixing 
depends only on perturbation parameter, consistent with 
experimental data. One may speculate that the back- 
ground and statistical errors may make these matrices 
random. 

In the following, we present the JPDF for the random 
reverse-cyclic matrices, and, an exactly solvable model 
related to this problem. The form of potential (for a 
single particle) has been discussed in quite a few phys- 
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ical situations. It has been interpreted as a screened, 
two-dimensional isotropic harmonic oscillator in a differ- 
ent context f22T |. It has found use in explaining roto- 
vibrational states in the case of diatomic molecules by 
conside ring a five-dimensional version of the Davidson os- 
cillator 1231 . and in a different context of dynamical sym- 
metries [24j and uncertainty relations [251 ] . In the con- 
text of many-body physics, we might imagine the above 
Hamiltonian describing bosons in a harmonic trap (r 2 
term), interacting via a dipolar electric field (r~ 2 term). 

We collect some known results related to the eigen- 
decomposition of reverse-cyclic matrices. A known eigen- 
decomposition becomes a very advantageous tool, to de- 
rive the joint probability distribution function for eigen- 
values. Karner et al. p6j have shown that the eigen- 
decomposition for an odd-dimensional reverse cyclic ma- 
trix is given by 



The JPDF in matrix space will be given by, using (J7J , by 
/o A \ (3/2) 

P(a,b,c) = exp[-3A(a 2 + 6 2 + c 2 )]. (9) 

From Eq. (J2), we can diagonalize H and it is also clear 
that there are only (n+l)/2 independent eigenvalues for 
odd-dimensional matrices. For the 3x3 case, the explicit 
form of R is 







R = 





^expH0/2) -*= exp (-i0/2) 
,(1 -±=ex P (*0/2) -^exp^/2), 



(10) 



A = (Ei, \E 2 \,..., \E(n_iy 2 \i -\E( n -l)/2\, 



(2) 



It takes a simple algebra then to show that 
0= -(Ei + 2|£ 2 |cos6>) 
b=-(Ei- \E 2 \ (cos8 + VZsmO 



(11) 



\E 2 \) 



c= - (Ei - \E 2 \ (cosO-Vs 



sin ( 



F rs (n) = —=e 



2iri(r-l)(s-l)/r, 



r,s = 1,2,..., n. (3) 



~ -iI k EI k ' 



(4) 



with = 



^diag f e^ 1 



/2 p i<Pk/2 



Ifc is an anti-diagonal identity matrix, and < (f>j < 27r. 
The eigen-decomposition for an even-dimensional reverse 
cyclic matrix takes following form, 

"-"(UK*)* < 5 » 

Al = (-El, |i?2|, • • • ) \E( n ^ 2 )/ 2 \, E n / 2 . -\E( n _ 2 y 2 \, ■ ■ ■ ) ~ l^l) 

with 

/£* iEt/ fc 
/,': 1 I e r-^ + 'x2/,- 



Using (TTTj) in ([5]), we can find the jpdf for eigenvalues 
and an independent parameter 9 coming from the eigen- 
vector. Note that in H, the independent parameters are 
three in number, namely a, b and c; while in the eigen- 
decomposition, we have Ei,E 2 ,9. The Jacobian for the 
transformation (fTTj) is given by 2 j^ . The jpdf for eigen- 
values is 

P(E X ,\E 2 \,9) = ^| (^j m) ex P [-A(El+2E*)] 

where E\ G (—00, 00), \E 2 \ € [0, 00), 
0€[O,27r). (12) 

Notice that the domain of \E 2 \ is [0,oo), and that the 
function on the right hand side is an even function of E 2 , 
thus we can rewrite the JPDF after an integration over 
9 in the following form, 



(6) 



P(Ei,E 2 ) = 2tt 



\E 2 



(-) 

3^3 v 7T y 



(3/2) 



II. JPDF, SPACING DISTRIBUTION AND 
DISCUSSION 

Consider an ensemble of reverse cyclic (RC) matrices, 
drawn from a Wishart distribution, 



exp(~A(E 2 1+ 2E 2 2 )). 

_ (13) 

The density of Ei (the trivial eigenvalue) [27( comes out 
to be Gaussian as expected because of Ei being a sum of 
Gaussians. On the other hand, the density of non-trivial 
eigenvalue E 2 is given by HI]). 



p(E) = 2,4|.E|exp(-2AE 2 ) 



(14) 



P(H) ~ exp (-ATr(H^H)). 



(7) 



Let us start with the simplest case, namely an ensemble 
of 3 x 3 reverse cyclic matrices, 



H 




(8) 



Also, due to product structure of the JPDF, the den- 
sity of non-trivial eigenvalue will remain the same for 
higher-dimensional matrices. The presence of \E\ ensures 
that there are no non-trivial eigenvalues present at origin 
while they increase linearly along both the positive and 
negative real axis. It is as if there is a hole in the density 
of non-trivial eigenvalues (see FigJT|. This has been inde- 
pendently derived by Bose et al. |15j without obtaining 
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FIG. 1. (Color online)Normalized density of non-trivial eigen- 
values for an ensemble of 20000 reverse cyclic matrices of size 
15 x 15 is compared with the analytical form. The density is 
normalized such that averaged density for positive eigenvalues 
is 1/2. 



the JPDF. Also notice that, it is the limiting distribution 
in the case of [l5[ while here it is an exact result for any 
dimension (matrix). The spacing distribution between 
Ei , E-x can now be calculated as 

/oo poo 
/ P(E 1 ,E 2 )5(s 12 -\E 1 -E 2 
-oo J —oo 

12VAe 



-oo 
-As? 



9^ 



AAe~ 2A ^ 



37rsi2Erf 



4*12 



9^ 



(15) 



The value of A can be chosen so that s\ 2 P(si 2 )dsi 2 — 
1. A numerical histogram is compared with (|15l) in FigJ5] 
One could think of spacing between the second and third 
eigenvalue of H, but due to their special form as \E 2 \ and 
— \E 2 \, it is simply given by s 23 = 2 j -E2 1 , so the spacing 
distribution as expected is very similar to the density of 
\E 2 \ and is given by (JT6J) : 



P{s2z) = As 23 e' 



2 s 23 



(16) 



Again, the value of A is chosen such that 
Jo s 2 3P(s 2 3)ds 2 3 = 1, which turns out to be ir/2. 
A comparison with the numerical data is shown in Fig. 
® 

In the case of 5 x 5, a similar procedure will give the 
JPDF as in dTTJ) with Ei £ (-00, 00) and 0, € [0, 2tt): 



P(Ei,E 2 ,E 3 ,6i,9 2 ) 



\Eo\\Ea 



exp 



25\/5 
A(El 



5A 

IT 



5/2 



2E 2 2 + 2E 2 3 )](17) 

The density of -EjS and the spacing distribution for the 
cases appearing in 3 x 3 reverse cyclic matrices remain 
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FIG. 2. (Color online) Spacing distribution S12 for an ensemble 
of 20000 reverse cyclic matrices of size 3x3, 5x5 and 15x15 
is compared with the analytical form (eq |15|) . 



the same. There is an additional spacing possible, namely 
between two positive \E 2 \ and I-E3I. Let us denote this 
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FIG. 3. (Color online) Spacing distribution S23 for an ensemble 
of 20000 reverse cyclic matrices of size 3x3, 5x5 and 15x15 
is compared with the analytical form (eq |16|l . 



by s pp . Its distribution is 

/oo />oo />oo 
/ / 4P(E 1 ,E 2 ,E 3 )S(s pp -\E 2 -E 3 \) 
-00 Jo Jo 

= As pp e~ 2As pp 

-±V^Ae- As ^{-l + 2As 2 pp )Eifc(VA Spp ). (18) 

The area under this distribution is 1/2. Taking care of 
the domains of \E 2 \ and I-E3I, and accounting for the 
spacing between these and between — \E 2 \ and — I-E3I, we 
obtain the correctly normalized distribution. This can 
be seen to be in agreement with the numerical data (see 
Fig. 2|). This same distribution [Eq. [THJhas been com- 
pared with the distribution of spacings among all posi- 
tive eigenvalues except the Gaussian distributed one of 
an ensemble of higher-dimensional reverse-cyclic matrices 
(e.g. 15 x 15). The agreement is good. For the (2n + 1) 
dimensional reverse-cyclic matrix JPDF is a straightfor- 
ward generalization of (fT7|) and is given by (p~9|) . 

/ A \ (2n+l)/2 

P{E\,E 2 , . . . , E 2n+ i,9i, . . . , 9 2n ) =1 — 1 

n+l 

\E 2 \... \E n+1 1 expl-AiEf + 2 ]T £?)](19) 

i=2 

This JPDF can be understood as follows. Clearly, the 
nature of the first (trivial) eigenvalue is very different 
from the others (nontrivial), and its distribution will be 
Gaussian. We focus on the rest of the eigenvalues. The 
diagonalizing equation, H = OAO T where O = R is 
an orthogonal matrix [261 ] , has the correct number of in- 
dependent parameters. For a (2n+l)-dimensional matrix 
H, the left hand side has only (2n + 1) independent vari- 
ables while the right-hand side has (n+l) independent 
eigenvalues with n angle variables in O. As dH will con- 
tain (2n+l) independent differentials, so a multiplication 
of H with a scalar a will satisfy d(aH) = a 2n+1 dH . Now, 
(n+l) of them will be absorbed in the scaling of measure 
dA[as independent eigenvalues are (n + 1)]. Hence, from 
the scaling property of d(aH), dH will be a homogeneous 
polynomial of degree n[28|. Our prototype examples for 
n = 3 and 5 has shown that they vanish linearly as eigen- 
values approach the origin, hence the polynomial in the 
eigenvalues is necessarily proportional to \E 2 \ . . . |.E n +i|. 
The even case is not very different from the odd one, ex- 
cept that E n / 2+1 appears along with E±, the rest being 
the same as that in (|19| . 

III. SCREENED HARMONIC OSCILLATOR 
AND JPDF 



Now we show the exactly solvable n-body problem, 
the ground-state wave function of which is such that 
the probability density has the same mathematical form 
as (fT9|) . It can be verified that IT9|) corresponds to 
|\&(xi, x 2l x n )\ 2 where ^(x\,x 2l ...,x n ) is the ground 
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FIG. 4. (Color online) Spacing distribution s pp for an ensem- 
ble of 5000 and 20000 reverse cyclic matrices of size 5x5 and 
15 x 15 is compared with the analytical form foq |18l) . 



state wave function with eigenvalue (4n — 3)^4 of the n- 
body problem with the Hamiltonian: 



U(xi,X2,—,X n ) 



a^ 1 + ±Ua^-^ 



i=2 



To illustrate that this is so, let us verify for n = 2. 
This will also be sufficient for general n due to the iden- 
tical form of separable H. As we need to take dou- 
ble derivatives of the wave function, it will be pru- 
dent to replace |_Ej|s by y/Ef- Hence, ^(xi,x 2 ) — 

c\Jy/x^exp {-A/2 [x\ + 2x1)), the JPDF for a corre- 
sponding three-dimensional reverse-cyclic matrix, 

d 2 



di 



rV(xi,X 2 ) = A(-1 



Axf) <&(x 1 ,x 2 ) 



-l + 16Ac| (-l + Axf)) 



cP_ 



A 2 x 2 + AA z x 



dx\ 



2„2 



^(x u x 2 ) 

l/(4z^ *{ Xl ,x 2 ) 
= 5A 1 $(xi,x 2 ) 



(20) 



This proves our assertion. 

This system has been the subject of a lot of work, 
initiated by Perelomov (29|. The only potential that can 
be added to a harmonic interaction is 2a /x 2 if we want 
to successfully construct the creation and annihilation 
operators for the above model [30| . This work relates 
this well-known model to a random matrix theory for 
reverse-cyclic matrices, which constitutes a remarkable 
addition to the known connections along similar lines. 

The linear level repulsion obtained here has its origin 
in the product of the absolute value of the eigenvalues 
in the JPDF. This is reflected in the interaction among 
eigenvalues if we write the JPDF as a partition function 
for an n-particle system. This interaction, in the con- 
text of random matrices is the Coulomb interaction in 
two dimensions. In contrast, the case of random-cyclic 
matrices has a JPDF which is just the exponential 
containing a sum of the square of the modulus of the 
complex eigenvalues. The eigenvalues are in a plane, and 
the level repulsion comes out as a Rayleigh distribution 
for the Poisson process on a plane, which has the same 
functional form as Wigner's spacing distribution for the 
orthogonal ensemble. Thus, we have a very interesting 
situation for the random reverse-cyclic and random cyclic 
matrices in that we obtain the same formula for the spac- 
ing distribution but the origin is different. 



IV. SUMMARY 

In summary, we have shown that reverse-cyclic ma- 
trices though a subset of symmetric matrices have an 
unusual density and spacing distribution. In contrast to 
semi-circle density, this ensemble admits a density with 
a hole at the origin. Again, the spacing distribution has 
a variety ranging from Gaussian-looking distributions to 
Wigner type distributions. We also observed that the 
JPDF is just the square of the modulus of the ground- 
state eigenfunction of an exactly solvable many-body 
Hamiltonian in one dimension, of a screened harmonic 
oscillator potential. Hence the correlations between the 
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different particles in the potential will be the same as sions. 
that derived from the joint probability distribution func- 
tion for the random matrix theory. 
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